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Abstract. There has been remarkable progress in numerical relativity recently. This has led 
to the generation of gravitational waveform signals covering what has been traditionally termed 
the three phases of the coalescence of a compact binary - the inspiral, merger and ringdown. In 
this paper, we examine the usefulness of inspiral only templates for both detection and parameter 
estimation of the full coalescence waveforms generated by numerical relativity simulations. To 
this end, we deploy as search templates waveforms based on the effective one-body waveforms 
terminated at the light-ring as well as standard post-Newtonian waveforms. We find that both of 
these are good for detection of signals. Parameter estimation is good at low masses, but degrades 
as the mass of the binary system increases. 
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1. Introduction 

Several ground-based interferometric detectors are now in operation to detect gravitational waves. 
These include the Laser Interferometric Gravitational-wave Observatory (LIGO) at two sites in 
Livingston and Hanford, USA, and the Virgo detector in Cascina, Italy. They have recently 
completed a first science run at or close to design sensitivity and are sensitive to gravitational 
waves from coalescing binaries at distances of tens to hundreds of megaparsecs depending on the 
total mass and the mass ratio of the system. The broadband sensitivity (40-400 Hz) of these 
detectors makes it possible to search for binaries with a rather large range of component masses 
from one to hundreds of solar masses. This range of masses includes both the neutron star binaries 
(which are known to exist) as well as neutron star-black hole and double black hole binaries (of 
which we have no observational evidence). 

In this article we test the efficiency of inspiral waveforms for the detection and parameter 
estimation of the full coalescence signal. We restrict our attention to two waveform families. 
The first of the families is the Fourier domain model, called TaylorF2 or SPA J^, which is an 
analytical approximation to the Fourier transform of the standard post-Newtonian (PN) [21 [3l IH [5] 
waveform (i.e., TaylorTS) computed using the stationary phase approximation. The highest PN 
order available in the LIGO Scientific Collaboration (LSC) code base [6 for this family is (f/c)^ 
(i.e. 3.5PN order). A number of searches by the LSC for compact binaries of low masses (i.e. 
M < 25 Mq) have used this model to second post Newtonian order as optimal templates [T, 'S^, '9', TO" . 

The second family we consider is the effective one- body (FOB) [UJ model at four P^S (i.e., 
(w/c)^) order but terminated at the light ring. As discussed below, the EOB model provides 
the full waveform, i.e., the inspiral, plunge, merger and ringdown waveform. In particular, the 
merger-ringdown waveform is obtained by stitching a superposition of quasi-normal modes to the 
inspiral-plunge waveform at the EOB light-ring [121 [T3l [TTl [Ml [TSl [T6] . In this study, however, we 
use the EOB model without the ringdown modesjj Thus, it captures some part of the coalescence 
signal, and is therefore better suited to search for higher mass signals. Here we will focus on the 
efficiency of our template bank to capture coalescence signals with the TaylorF2 and EOB models. 
In Section |21 we discuss in greater detail the dynamics of binary black hole mergers and the PN 
and EOB models. 

We test the efficiency and parameter estimation accuracy of the searches in two different 
ways. First, in Section [3] we perform a Monte-Carlo study of the efficiency of the TaylorF2 and 
FOB models to detect the full waveforms. Since the full waveform is not known over the entire 
parameter space, we make use of the EOB waveform calibrated to the numerical relativity results 
[14[ I15[ I16[ I17[ I18j . Then, in Section [4l we perform a similar comparison making use of waveforms 
generated numerically. This study was performed on the NINJA [19 data set which comprised 
simulated data for the LIGO and Virgo detectors with numerically obtained binary coalescence 
signals added. The original numerical results for the NINJA numerical waveform contributions are 
described in[2ni[2Il[ll[13[21[13[M[13[lHl[Il[Ml[3ni[3Il[3a (where these are pubhshed results) , 
the codes are described in[Ml[Ml[Ml[13[M[371|3Hl[Mllini[lHlllIllll[ll[32]- 

I Post-Newtonian expansion is currently known only to 3.5 PN order [21 |3l U [5] ; the unknown 4 PN coefficient is 
determined by matching the EOB to numerical relativity waveforms 1121 1131 1111 1141 1151 1161 (see below) . 
§ Throughout this article, we will use EOB to refer to the EOB model truncated at light ring, and EOBNR to refer 
to the EOB waveform including the ringdown modes, calibrated to numerical relativity results. 
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2. Binary black hole dynamics 

The evolution of a black hole binary is driven by back-reaction due to the emission of gravitational 
waves which causes the system to inspiral and merge. Confident detection of the emitted signal is 
greatly facilitated by an accurate understanding of the dynamics of the binary and the shape of the 
emitted waveform during inspiral and merger. The early evolution of a binary can be well-modelled 
by the PN approximation during which the system slowly inspirals on an adiabatic sequence of 
quasi-circular orbits located at the (stable) minimum of the changing effective potential. In fact, 
for most of its lifetime a binary black hole can be accurately described by the balance of the rate of 
change of the binding energy with the energy carried away to infinity by the radiation as given by 
the quadrupole formula. In fact, many of the ideas (effective potential, last stable orbit, etc.) that 
are relevant when the component masses are greatly separated are still very useful in analytically 
modelling the system close to coalescence. However, they are perhaps not so useful or needed from 
the view point of numerical evolution. 

As the system evolves, the effective potential changes and reaches a point when the potential 
transforms from one having a stable minimum and an unstable maximum to one having just an 
unstable minimum. After this, the system no longer possesses any bound orbits. The transition 
point, called the last stable or biiO (LSO), occurs when the radius r of the orbit (in Schwarzschild 
coordinates) approaches r ~ 6GM/c^, where M is the total mass of the system. In terms of the 
dominant component of the emitted radiation, this corresponds to a gravitational-wave frequency 
of /insp — 44OHz(A//lOM0)~^. Therefore, for masses less than about 10 AIq, only the inspiral stage 
of the coalescence lies in the detector's sensitive band of 40-400 Hz. 

Once the system passes the LSO, the two black holes plunge towards each other and merge 
in about one orbital time scale of the LSO to form one single distorted black hole. This is the 
so-called merger phase which is amenable to analytic description by a clever re-summation of 
the PN approximation but more recently numerical relativity simulations have provided a better 
understanding of the merger phase and continue to provide new insights. The frequency of the waves 
during this phase changes rapidly from /merge ^ 440Hz(M/10 AfQ)-^ to 1200 Hz (A//10 Afo)"^ 
During the late stages of the coalescence, the highly distorted black hole, that results from the 
merger of the two parent black holes, settles down to an axi-symmetric quiescent state by emitting 
its deformation in the form of gravitational waves. The radiation from this phase is well described 
by black hole perturbation theory and consists of a set of quasi-normal modes (often referred to as 
ringdown signal) whose fundamental frequency is /ring ^ 1800 (Af/10 A/Zq)"^ when two equal mass 
non-spinning black holes merge to form a single black hole whose spin magnitude is estimated to 
be J/M'^ ~ 0.7. The first two overtones of this mode have frequencies of ^ 1650 Hz and ^ 1700 Hz, 
for the same system. 

2.1. Search templates 

The foregoing discussion hints that binaries whose total mass is less than about 10 Mq can be 
detected by using templates that are described by the PN approximation. In fact, experience 
suggests that we could make do with the PN waveform as templates even when the total mass is 
as large as about 25 Mq and they have been used in the search for low-mass systems (i.e. systems 

II The LSO discussed in this paragraph refers to the test-particle hmit of a binary when the mass ratio is very small; 
PN corrections and resummed models modify the location of the LSO ^44 . However, these changes are unimportant 
to the present discussion where we are concerned with approximate numbers to determine the rough boundary 
between where a certain phase of the evolution is dominant. 
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with their total mass less than about 25A/0) in the data from LIGO and Virgo [45^. However, for 
higher mass black hole binaries (i.e. binaries with their total mass greater than about 25Mq) the 
merger of the binary occurs in the detector's sensitive band. At merger, the dynamics is no longer 
adiabatic and is, therefore, not well-modelled by PN expansion. It has been a long standing aim 
of numerical relativity to generate the full waveforms for gravitational wave detection from higher 
mass black holes. 

There has been significant progress recently in numerical relativity with several groups having 
successfully simulated the merger of two black holes (see, for example, [42l [39l [46] , and the NINJA 
related numerical relativity results cited earlier). The longest of these simulations last for tens of 
orbits |30j , and they cover different mass ratios and are beginning to explore the space of component 
masses with spins. 

Nearly a decade ago, analytical work by Buonanno and Damour extended the PN dynamics 
beyond the last stable orbit to calculate the merger dynamics. This analytical method, called 
EOB computes the dynamics up to the light ring of the effective potential and the waveform can 
be computed for separations larger than about r ~ 2.2 M. In this work we have used the EOB 
waveform terminated at the light ring as search templates. The EOB formalism, however, provides 
the full waveform, including the merger-ringdown portion which is attached at the EOB light ring. 
Moreover, the availability of numerical relativity simulations has helped in fixing certain unknown 
higher order (4PN) terms in the EOB model by fitting the analytical waveform to numerical 
relativity. The current implementation of the full EOB model in LAL uses the 4PN order parameter 
obtained in [17] by calibrating the EOB model to the Goddard-NASA numerical simulations, and 
three ringdown modes (i.e., the fundamental mode and two overtones). We note that, more recently, 
the EOB model has been further improved by calibrating it to more and more accurate numerical 
simulations [Ml [151 [El \lM which will form the basis of future searches. In the current study, 
EOBNR is used to calibrate the efficiency of the inspiral models. 

Ajith et al [47l have used a phenomenological approach to match the inspiral phase from 
PN approximation to the merger and ringdown from numerical relativity. Here, an hybrid waveform 
obtained by stiching together the PN inspiral waveform to the numerical-relativity inspiral, merger 
and ringdown waveform is first built. Then, the phenomenological waveform in the frequency 
domain is constructed by requiring closeness to the hybrid waveform. In the long run, it is likely 
that these full waveforms will be used as templates to search for inspiral signals in gravitational 
wave detectors. 

In order to test the efficiency of EOB and TaylorF2 families, we will use the EOBNR waveform 
as our "true" waveform and see how well these partial waveforms perform in both detection and 
parameter estimation. 

2.2. An example waveform 

Fig. [T] shows the waveform expected from a pair of non-spinning 10 Mq black holes during the last 
50 ms before merger. The left panel shows the time-domain waveform h(t) and the right panel 
compares the time-domain waveform with the signal 'perceived' by the initial LIGO detector whose 
noise has been whitened. In other words, the right panel plots whitened template q{t) given by 
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Figure 1. An example of the EOBNR waveform for a binary consisting of two equal mass, non- 
spinning black holes each 10 A/q . The EOB dynamics allows the computation of the inspiral (left 
panel, black solid line) and plunge (left panel, roughly the first two cycles of the red dashed line) 
phases but the merger and ringdown waveform (left panel, latter part of the red dashed line) is 
stitched to the end of the plunge phase by matching the amplitude of the waveform and its first 
two derivatives by using a superposition of the fundamental quasi-normal mode and its first two 
overtones. The right panel compares the time-domain signal (red dashed line) with the whitened 
signal (blue solid line). 



where H{f) — h{t) exp(27ri/i) At is the Fourier transform of the time-domain signal and Sh{f) 
is the one-sided noise power spectral density of initial LIGO. In the left panel, the inspiral part of 
the signal given by EOB dynamics is shown in black solid line followed by the plunge, merger and 
ringdown phases in red dashed line. The right panel compares the time-domain waveform in the 
left panel (red dashed line) with the whitened signal (blue solid hne). 

Note that although the time-domain signal is dominated by the plunge, merger and the 
ringdown phases, the detector noise spectral density (i.e, Shif)) suppresses them, making the 
inspiral phase more dominant. For systems with greater masses, more of the merger phase appears 
in band. For systems with total mass larger than about 80 Mq the merger and ringdown signals 
begin to dominate over the inspiral phase. For such systems it is important to deploy EOBNR 
templates. As a result, we cannot expect our template families to do well in capturing high mass 
binaries. 

3. Bank Efficiency 

Matched filtering, the data analysis technique used in most searches for binary black holes, is pretty 
sensitive to the phasing of the signal, which in turn depends on (a subset of) the source parameters. 
In the case of non-spinning black holes on a quasi-circular orbit, the only parameters that we must 
consider are the two masses. The location of the binary on the sky, the distance to the source, 
the polarization of the wave, etc., are not important for a single detector as they simply affect the 
amplitude of the waveform. Although we won't know the time at which the binary merges nor the 
phase of the signal at that epoch, these parameters need not be explicitly searched for [I], and are 
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easily extracted in the process of maximising the cross-correlation of the template with the data. 



3.1. Template bank 

Our goal in this Section is to study the efficiency of the two template families in detecting binary 
black hole coalescences. To this end we first set up a template bank — a set of points in the 
parameter space of the component masses. A geometric algorithm described by Babak et al 
is used to generate the template bank and it is the same algorithm irrespective of which family 
of waveforms is used to filter the signals. The bank is designed to cover the desired range of 
component masses of the binary. In addition to the range of the component masses, our template 
bank algorithm requires us to specify a parameter called the minimal match, MM. The minimal 
match is the smallest overlap guaranteed between a signal with random source parameters and the 
template nearest (in the geometrical sense) to it in the parameter space. 

Assumptions made in the construction of the template placement algorithm imply that this 
will be strictly true only when (a) the templates and signals belong to the same family, and (b) 
the ending frequency (i.e., the LSO or light ring depending on the waveform in question) is greater 
than the upper end of the sensitivity band. The latter condition further implies that we can hope 
to achieve overlaps of MM or greater only for waveforms whose total mass is smaller than a certain 
value depending on the detector bandwidth; in the case of initial LIGO this is 10 Mq. We have 
chosen MM = 0.95. However, since our templates and signals belong to different PN approximations 
we cannot expect to achieve this overlap even for total mass less than 10 Mq. 

The template placement algorithm chooses a hexagonal grid in the two-dimensional parameter 
space of the component masses and it is an optimal algorithm in the sense that it gives the smallest 
number of templates possible for a given minimal match [50j . 

Having constructed a template bank we then generate an EOBNR signal with the values of its 
masses, the epoch and the phase at coalescence, all chosen randomly but in a given range. Since 
the EOBNR extends beyond the LSO, it is possible to generate signals with the total mass in the 
range [10, 300] and minimum component mass of 5Mq. We performed simulations for both 
the TaylorF2 and EOB templates over a restricted range of the parameter space and an additional 
EOB analysis over a much larger range. Specifically, in the case of TaylorF2 we set up a template 
bank with the total mass in the range [6, 80] Mq, with minimum component mass of 3Mq, and 
searched for signals with total mass in the range [10, 72] Mq. Although the TaylorF2 model is 
not expected to be a good approximation above about ~ 25 Mq, we are using it to much higher 
mass. For the EOB family truncated at light ring, we performed two separate analyses. The first 
had templates covering the range [6, 120] Mq, minimum component mass 3Mq and injected signals 
with total mass in the range [20, 80] Mq. Since the EOB waveforms extend to a significantly higher 
frequency than the standard post Newtonian ones, we also performed an analysis with templates 
of total mass in the range [60, 400] Mq and minimum component mass of 30Mq. In this case the 
template bank consisted of only 9 templates. The simulated signals had total mass between 70 and 
300Mq. 

Next, for each point in the template bank we generate waveforms from our template families 
(EOB and TaylorF2) and measure their overlap with the random signal. The overlap Ok of the fcth 
template qk{t; rui, mf) and the signal h{t) is defined by 

/" , , , , _ . .. d/ 



Okimt ml) = max2 / [H{f)Ql{f- m\, m%) + H*{f)Qk{f; m^)] e 



-27ri/i 
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where H{f) and Q{f) are the Fourier transforms of h{t) and q{t), respectively, and Q* is the 
complex conjugate of Q. This allows us to compute the maximum overlap between our template 
waveforms and a random signal^. This process is repeated for 1,000 different realizations of the 
random mass parameters and the maximum of the overlap over the entire template bank is recorded 
in each case. We will now discuss the results of these simulations. 

3.2. Efficiency for detection 



EOBNR injections EOBNR injections 
I.n r- ' ' ' ' ' ' 1 1.00 I • ' • 




Injected Total Mass {Mq) Injected Total Mass (Mq) 

Figure 2. These plots depict the efficiency of the TaylorF2 (left panel) and BOB (right 
panel) template banks in detecting the coalescence waveforms assumed to be well-represented 
by EOBNR. Each dot corresponds to the overlap of a random EOBNR signal maximized over the 
template bank consisting of TaylorF2 or EOB waveforms. 

Fig. [2] plots the results of our simulation. The left panel shows the (maximum) overlaps of the 
TaylorF2 template bank with random signals, one dot for each trial. The right panel shows the 
same but for the EOB template bank. The TaylorF2 model has overlaps larger than 90% for only 
signals whose total mass is less than about 22 Mq. The overlap falls off quickly for masses larger 
than this, reaching slightly more than 0.5 when the total mass is about 50 Mq. The fact that the 
overlap remains unchanged beyond 50 Mq is probably due to the TaylorF2 template matching the 
later merger or ringdowm part of the EOBNR waveform. In contrast, the EOB model achieves 
an overlap greater than 90% over the entire mass range. This is due to the fact that the EOB 
waveform extends to the light ring. 

Figure [3] shows the overlap obtained using the high mass EOB template bank and an extended 
range of EOBNR simulated waveforms. We obtain overlaps of better than 85% for systems whose 
total mass is less than 150 Mq, and the overlaps remain more than 65% even for systems with total 
mass less than 200 M©. The oscillatory behaviour seen in this case is due to an edge effect: the 
density of templates gets smaller and smaller as the total mass becomes larger. For this range of the 
total mass, there are only a handful of templates in the entire bank; as a result, a single template 
might be available for a pretty large range of masses, causing the overlaps to swing up and down 
as the total mass is increased. 

^ For the sake of saving space we have not discussed the maximization over the phase of the signal. This can be 
found, for instance, in Ref. [l]. 
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Figure 3. This plots depict the efficiency of the EOB template bank generated for the mass 
range [60, 400] Mq in detecting the coalescence waveforms assumed to be well-represented by 
EOBNR. Each dot corresponds to the overlap of a random EOBNR signal maximized over the 
template bank consisting EOB waveforms. 



The overlaps are surprisingly high, especially considering that there are only nine templates 
in the bank. It is possible that the very high mass EOB signals are very short in the detector 
band (perhaps a cycle or two) and the abrupt cutoff of the template as a result of termination at 
the light ring could bleed power into a frequency region where there is no power in realitjEI- This 
spurious template power may lead to large overlaps with the merger and ringdown parts of the 
EOBNR waveforms. TaylorF2 does not suffer from this predicament. This is because TaylorF2 is 
generated in the Fourier domain and the abrupt cutoff of the signal does not cause any problem 
in the frequency domain and we are unconcerned with spurious effects in the time-domian as they 
occur outside the region of our interest. 

3.3. Efficiency for parameter estimation 

The matched filtering statistic gives the likelihood for a signal to be present in the data as opposed 
to the data being pure background noise. The parameters of the template which maximize the 
likelihood are maximum likelihood estimates. Having determined the efficiency of our template 
banks in capturing the coalescence signals, we next consider how good they are in measuring the 
signal parameters in the maximum likelihood sense. If the template waveforms and the signal they 
are intended to detect both belong to the same family then in the limit of large signal-to-noise ratio 
the distribution of the maximum likelihood estimates will be centred on the true signal parameters. 
We are in a situation wherein the template waveforms and the signals they are intended to recover 
belong to different approximations. Therefore, one can expect a systematic bias in the estimation 
of parameters. 

To gauge the reliability of the two families in estimating the parameters of the true signal we 
make use of the results of the simulation from the previous Section. This simulation computed the 

+ The spurious power is in itself not a bad thing but large noise glitches in the region where there is spurious power 
could cause false alarms. This is especially the case when the detector noise is contaminated by large amplitude 
non-stationary noise glitches. 
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overlap of the templates with the signals in the absence of any noise. Therefore, the parameters of 
templates that maximized the overlap when compared to the true parameters of the EOBNR signal 
give a measure of the systematic bias in parameter estimation due to the difference in the waveform 
families representing the templates and the signal. Fig. |4] shows the measured total mass of the 
template (y-axis) vs the true total mass of the EOBNR signal (x-axis) and the colour is determined 
by the corresponding overlap. 
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Figure 4. Parameter estimation accuracies for tlie TaylorF2 (left panel) and the EOB (right 
panel) models. The total mass of the injected EOBNR signal is plotted vs the total mass of 
the template that obtained the best overlap. The colour of the plotted point is determined by 
the overlap between signal and template. Clearly, there is a positive correlation between the 
injected and measured masses in the case of TaylorF2 model, with the spread in the measured 
values becoming larger at higher masses. Additionally, for higher masses the overlap decreases. 
Interestingly, all simulated signals with the total mass greater than 50Mq are recovered with the 
same template of mass 35Mq . The EOB templates show a similar correlation between simulated 
and recovered masses. However, since the EOB templates extend to higher frequency, the EOB 
templates are capable of accurately recovering the higher mass simulations with accurate mass 
parameters and an overlap greater than 90%. 



Both the TaylorF2 and EOB models provide good parameter recovery as well as good overlap 
for simulated signals with a total mass less than 30Mq. At masses higher than this, the TaylorF2 
waveforms achieve a significantly lower overlap, although parameter estimation remains good up to 
a total mass of about 50Mq. Above this mass, all simulated signals are recovered by a template of 
mass 35Mq. Interestingly, this template will terminate around 120Hz which is close to the most 
sensitive frequency of the LIGO detector. Thus, it seems likely that the inspiral template is picking 
up power from the merger and ringdown parts of the EOBNR waveform. For the EOB templates, the 
overlap does drop off somewhat for higher masses, but the parameter recovery remains reasonable 
throughout, with a slight tendency to underestimate the total mass of the signal. For higher masses. 
Figure [5] shows the parameter estimation accuracy for the EOB waveforms. Despite the fact that 
there are only a handful of templates in the bank, the parameter recovery is reasonable. 
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Figure 5. Parameter estimation accuracies for high masses using the EOB templates. The total 
mass of the injected EOBNR signal is plotted vs the total mass of the template that obtained 
the best overlap. The colour of the plotted point is determined by the overlap between signal 
and template. Even for high masses, there is a good correlation between the simulated and 
recovered masses. The discreteness of the EOB template bank is clearly seen from the limited set 
of recovered masses, and indeed this seems to account for the majority of poor overlaps. 



4. NINJA Results 

The Numerical INJection Analysis (NINJA) project was a mock data challenge, where the data were 
generated at the design sensitivity of the initial LIGO and Virgo detectors and numerical relativity 
waveforms provided by a number of groups were added to the data. A number of data analysis 
methods were applied to the data, and the results of the NINJA project are available elsewhere 
[T^. For the NINJA analysis, we performed several runs through the data using the LSC's Compact 
Binary Coalescence (CBC) analysis pipeline. Here, we restrict our attention to two runs through 
the data which are similar to the TaylorF2 and EOB analyses described in the previous Section. 
This allows us to investigate the issues of detection and parameter estimation using these templates 
to search for full waveforms obtained from numerical relativity. The results are similar to those 
obtained in the previous Section, namely that the EOB search has a greater efhciency than the 
TaylorF2 search, but that both can detect high mass signals, although the parameter estimation is 
poor. 

The CBC pipeline was designed to analyze data from a network of detectors to search 
specifically for gravitational wave signals from binary neutron stars and black holes [9J . It proceeds 
as follows: First, a bank of templates covering the desired mass range is produced. For the NINJA 
analysis we used a template bank covering masses between 20 and 90Mq with a minimal match 
of 0.97. The data from each of the detectors is separately match filtered against the template 
waveforms [5T, . For this analysis, we restricted attention to the mock data generated for the three 
LIGO detectors (the 4 km and 2 km detectors denoted HI and H2, respectively, at Hanford, WA 
and the 4 km LI detector at Livingston, LA). A trigger is produced whenever the signal-to- noise 
ratio exceeds the desired threshold of 5.5. A coincident trigger is recorded whenever there are 
triggers from two or more detectors with comparable masses and coalescence times |52| . Finally, 
these coincident triggers are subjected to a set of signal based vetoes, in particular the [53] and 
[54j tests, designed to separate signals from non-stationary transients in the noise. 
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Template 


TaylorF2 


EOB 


Freq. Cutoff 


LSO 


Light Ring 


PN Order 


2 pN 


2 pN 


Found Inj, Single Detector (HI, H2, LI) 


72, 43, 66 


91, 64, 82 


Found Inj, Coincidence 


59 


83 


Found Inj, Coincidence + Signal Vetoes 


59 


80 



Table 1. Results of inspiral search for NINJA waveforms. There were 126 injections performed 
into the data. The table above shows the number of injections which were recovered using the two 
waveform families. The EOB search shows a significantly higher sensitivity than the TaylorF2 
waveforms evolved to LSO. Note that virtually all simulations which pass the initial coincidence 
requirement also survived the signal consistency checks. 
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Figure 6. Found and missed numerical injections for the TaylorF2 (left panel) and EOB (right 
panel) searches of the NINJA data. The found and missed injections are plotted on the total 
mass, Hanford effective distance plane. The effective distance for a detector provides a measure 
of the amplitude of the signal at that site, taking into account the distance and orientation of the 
source. For both searches, the majority of the close simulations are recovered. EOB templates 
are seen to perform better, particularly at higher masses. In both searches, the simulation with 
total mass 105Mq and effective distance of 200Mpc is missed. This has a peak amplitude at the 
start of the waveform, rather than coalescence, and although a coincident trigger is recorded the 
time between it and the waveform peak is too large for our algorithm to associate them. 



Our results are summarized in Table [T] which shows the number of injections recovered by the 
analysis pipeline at each stage of the analysis for the two searches described here. The EOB search 
is capable of detecting a greater fraction of the simulated signals than the TaylorF2 templates 
truncated at LSO. This is further highlighted in Figure [S] where we show those simulated signals 
which were recovered by the two different waveform families. The EOB model clearly performs 
better, particularly at higher masses. This is consistent with the findings of the previous Section. 

Next, we turn to parameter estimation. Figure [7] shows the accuracy with which the total 
mass of the simulated signals is recovered using the inspiral only waveforms. For both the TaylorF2 
and EOB models, the parameter recovery is poor, particularly at higher masses. This is to be 
expected, since we are searching with partial waveforms and, at the higher masses, it is the merger 
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Figure 7. Accuracy of recovering the total mass of simulated signals for the TaylorF2 (red +) 
and BOB (blue x) models. For both of the searches, the total mass is estimated poorly and 
systematically lower than the simulated mass. This is due to the fact that the search has been 
performed with inspiral only waveforms for which the search extends only up to 90Mq . 



and ringdown of the simulations which occupies the sensitive band of the detectors. Furthermore, 
the template bank extends only to a total mass of 90AfQ making accurate parameter recovery of 
the high mass simulations impossible. 

5. Discussion 

In the coming years, the first detection of gravitational waves from coalescing binaries will surely 
be achieved. Following the first detection, attention will focus on extracting as much astrophysical 
information as possible from the observed signal. In the studies described here, we have addressed 
the ability to perform both the detection and parameter estimation problems using template 
waveforms which cover only part of the binary coalescence. We have made use of two different 
waveforms — the TaylorF2 PN waveforms taken to LSO and the FOB waveforms to light ring. In 
addition, we have used two different methods to evaluate the detection and parameter estimation 
capabilities of these signals — a Monte-Carlo study using FOBNR waveforms as the "true" signal, 
and an analysis of numerical relativity waveforms in the NINJA data. In all cases, the conclusion 
is the same: the inspiral only templates are useful for detection of the signal, but do not provide 
good parameter estimation, particularly for the higher mass signals. This is to be expected as, for 
high mass binaries, it is the merger and/or ringdown which occurs at the most sensitive frequency 
of the detectors. We observe that the FOB waveforms perform somewhat better than the TaylorF2 
waveforms. This is expected as the FOB waveforms extend to the light ring and therefore capture 
the plunge part of the waveform which is not incorporated in the TaylorF2 PN model. 

The results of this study show that post-Newtonian based inspiral only waveforms will not be 
sufficient for satisfactory detection and parameter estimation of higher mass black hole binaries. Full 
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waveforms derived from a synthesis of post-Newtonian waveforms and numerical relativity results, 
such as the EOBNR model [Ml [15l [HI [ITl [T8j , or phenomenological models, such as [43|48], will 
be necessary. 
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